Dynamics of electric transport in interacting Weyl semimetals. 
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The response to an electric field (DC and AC) of electronic systems in which the Fermi "surface" 
consists of a number of 3D Weyl points (such as some pyrochlore iridates) exhibits a peculiar combi- 
nation of characteristics usually associated with insulating and conducting behaviour. Generically a 
neutral plasma in clean materials can be described by a tight binding model with a strong spin-orbit 
interaction. A system of that type has a vanishing DC conductivity; however the current response 
to the DC field is very slow: the current decays with time in a powerwise manner, different from an 
insulator. The AC conductivity, in addition to a finite real part a' (Q) which is linear in frequency, 
exhibits an imaginary part a" (Q.) that increases logarithmically as function of the UV cutoff (atomic 
scale). This leads to substantial dielectric response like a large dielectric constant at low frequencies. 
This is in contrast to a 2D Weyl semimetal like graphene at neutrality point where the AC conduc- 
tivity is purely pseudo-dissipative. The Coulomb interaction between electrons is long range and 
sufficiently strong to make a significant impact on transport. The interaction contribution to the AC 
conductivity is calculated within the tight binding model. The result for the real part expressed via 
the renormalized (at frequency Q) Fermi velocity v is: Act' (Q) — e*Q/ (9Tv^hv) [2 log (fi/Sl) — 5]. 

PACS numbers: 72.80.Vp, 11.30.Rd, 11. 15. Ha 



INTRODUCTION AND SUMMARY 

A long time ago a rather unorthodox physics of crys- 
tals possessing three-dimensional (3D) pseudo - relativis- 
tic quasiparticles [l|, exhibiting an electronic dispersion 
relation £k = w |k|, where the velocity v is of the order of 
the Fermi velocity in regular condensed matter systems, 
was invoked to describe properties of Bi. The ultra- 
relativistic linear dispersion relation describes two coni- 
cal bands (of opposite orientation) sharing the same cone 
tip. Recently several proposals 12tJ: \% ISf revived an in- 
terest in materials with such excitations nowadays called 
Weyl semi - metals. The Fermi "surface" of such materi- 
als, typically with dominant spin-orbit interactions, con- 
sists just of a finite number of disconnected points (called 
Weyl or Dirac points, defined below) rather than forming 
a continuous Fermi surface like electrons in usual metals. 
The revived interest emerged of course after years of in- 
tense experimental and theoretical study of graphene, a 
2D Weyl quasiparticle material [7|. Suspended graphene 
is just such a "semimetal" system and exhibits a number 
of remarkable properties. For example, despite having 
zero density of states at the Fermi level and ideally no 
impurities, it still has a nonzero DC conductivity |8|. 

While this band " touching" at a singular point was no- 
ticed in band structure calculations even before the semi- 
nal work of Wallace on graphitel9J, the first explicit use of 
the Dirac model in 3D was in the context of the two-band 
approximation model of bismuth[l|. Due to the strong 
spin-orbit interaction, the linear in k terms in a low en- 
ergy effective theory near the crystallographic point L of 
the FCC Brillouin zone, dominate over quadratic terms 
(that dominate near the F point, leading to a common 
"effective mass" description). The electronic excitations 



are described by an analog of the Weyl equation of parti- 
cle physics, which describes eight two-component chiral 
spinors "0 (two for each of the four L points) 



dttp± = ±wcr • VV'i. 



(1) 



The sign describes the chirality of the mode. Metallic 
bismuth is only approximately described by the ultrarel- 
ativistic "massless" dispersion relation since the quasi- 
particles of the opposite chirality are coupled and form 
four-component massive Dirac bispinors. In Bi therefore 
electrons are not Weyl, but lOmeV massive Dirac elec- 
trons where, in addition, the Fermi level is located away 
from the Dirac point. 

A number of related suggestions for suitable realiza- 
tions of Weyl semimetals were recently put forward. 
Kariyado and Ogata [SJ calculated the band structure of 
cubic inverse perovskites like Ca^PbO with significant 
spin-orbit coupling. They observed the appearance of 
six Weyl points with a very small relativistic electron 
mass down to AmeV on the line connecting the F - and 
the X - points. In iridium-based pyrochlores such as 
Y2lr20T, there arc Nw = 24 Weyl points located near 
the four L points of the FCC lattice. As noted in[2| these 
materials "in particular provide a unique opportunity to 
study the interplay of Coulomb interactions, spin-orbit 
coupling, and the band topology of solids" . Also, strong 
spin-orbit interactions can lead to a novel phase of mat- 
ter, the topological insulator[5| and various possibilities 
to create Weyl fermions combined into coincident oppo- 
site chiralitypoints or separated in the Brillouin zone in 
Bi02/Si02m, NasBi, HgiCr2Se4^. These proposals 
generated a great deal of experimental efforts [10|. The 
system with 3D Weyl points was proposed to appear in 



optical lattices [ll| following the discovery of "artificial 
graphene" [12| . 

Since the density of carriers in 3D Weyl semimetals at 
zero temperature is zero (as in suspended graphene in 
2D), the Coulomb interactions are unscreened and there- 
fore are expected to be important to the understanding 
of the electrical and optical response of these materials. 
Unsophisticatcdly the dimensionlcss coupling, 
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(2) 



is of order 1, provided the dielectric constant e is not 
large, since the analog of the light velocity, v, is of the 
order typically of the Fermi velocity. Note that the same 
Coulomb potential l/r created by an electron influences 
many more electrons in 3D compared to 2D, so naively, in 
3D its importance is expected to increase. While electric 
transport in noninteracting 3D Weyl fermions was stud- 
ied ( ISJ, [ij] and references therein to earlier papers in 
the context of particle physics), the contributions of po- 
tentially very important Coulomb interactions (Coulomb 
scattering corrections to transport as an example) have 
not been studied theoretically, except basic renormaliza- 
tion effects 1J|. This is in contrast to the situation in 
graphene. 

The effect of the Coulomb interactions in undoped 
graphene turned out to be highly nontrivial, even 
within perturbation theory, and have evoked a scien- 
tific controversy [iSfjlSj due to the problem of the "ultra- 
violet regularization" , and was just recently settled[19| 
by noting that the ambiguities are associated with the 
treatment of the separation of scales related to the chi- 
ral anomaly [20J. Some aspects of the Weyl semimetal 
physics are not dominated by the Weyl points of the 
Brillouin zone at which the spectrum is gapless. For ex- 
ample, the AC conductivity of undoped graphene (the 
weak logarithmic renormalization of the electron veloc- 
ity [2l| does not influence the result), is given in terms 
of its value in the noninteracting theory, ag = e^/4?i, by 






= l + Ca + (a^) 



(3) 



and is strictly dissipative (real). This expression is in- 
dependent of frequency (provided corrections of order 
hfl/j, 7 = 2.7ey being the hopping energy, are ne- 
glected). The value of the only numerical constant C 
appearing here has been a matter of intense controversy. 
The first detailed calculation [15| utilizing a sharp mo- 
mentum cutoff regularization of the Dirac model pro- 
vided a value C^^) = ft - f «== 0.51 of order 1. The 
use of the sharp momentum cutoff was criticized by 
Mishchenko[16J, who obtained an exceptionally small 
value of Ct^) = i| _ | « O.Ol making a "soft" mo- 
mentum cutoff regularization. He supported this choice 
by the consistency of the Kubo and the kinetic equation 



calculations of conductivity with that of the polarization 
function. The consistency required a modification of the 
long-range interaction so that it becomes UV cutoff de- 
pendent. This apparently closed the issue. Albeit such a 
small numerical value would have profound physical con- 
sequences even beyond the transport and dielectric prop- 
erties. Nevertheless the interaction strength C was recal- 
culated once again by Vafek, Juricic and Herbut[17[, who 
argued that the modification of the interaction requires 
simultaneously a Pauli-Villars regularization of masslcss 
fermions. They applied yet another regularization, mak- 
ing the space dimensionality fractional, D = 2 — e, that 
modified both the current operator and the interaction in 
such a way that they satisfy the Ward identities and ob- 
tained C^'^-' ~ T~f ^ 0.26. The dimensional regulariza- 
tion is questionable on physical grounds and in a compre- 
hensive subsequent work|18| the authors reaffirmed the 
small value C^^-*; it seems that this value is the commonly 
accepted one. The tight binding calculation [19[ however 
demonstrated that C*-"^^ is the correct one. To reveal 
the origin of the ambiguity exhibited by the various val- 
ues of C, the authors made use of a dynamical approach 
developed earlier [2^ (to address the problem of separat- 
ing the interband contributions from the intraband ef- 
fects due to contacts [23|,l2J|) directly in the DC case by 
" switching on" a uniform electric field in the tight binding 
model with Coulomb interactions, and then considering 
the large-time limit. This approach (known in field the- 
ory as the " infinite hotel story^M is the best way to reveal 
physical effects of anomalies j25|. One can directly sep- 
arate the contributions from the neighborhood of Dirac 
points and the "anomalous" contributions from the rest 
of the Brillouin zone, so that one can decide what regu- 
larization of the effective Weyl theory is the correct one. 
In this sense this is advantageous over the standard di- 
agrammatic Kubo formula calculation in continuum [l5|- 
|l8[ that might encounter the so-called Schwinger terms 
(found in Quantum Electrodynamics, that is similar to 
the 3D Weyl fermions). 

The purpose of the present paper is to study the effect 
of Coulomb interactions in 3D Weyl fermion systems with 
emphasis on dynamical aspects of electric transport and 
compare/contrast it with the corresponding results in 2D 
Weyl fermions. To achieve this goal we define in section II 
a tight binding description of Weyl fermions on a hyper- 
cubic lattice of any dimensionality similar to a variant of 
the Hamiltonian lattice model in field theory (25J alr eady 
used in its Lagrangian version to simulate graphene [26 1. 
In this model the electron's spin is strongly coupled to 
momentum and therefore the model is very reminiscent, 
in this respect, of the Wolff model of bismuth or lattice re- 
alizations of topological insulators. The lattice "regular- 
ization" is necessary to address the ultraviolet divergen- 
cies at the intermediated stages of calculations, a problem 
mentioned above. The universality of this description of 
Weyl fermions (of various origins) is supported by the 



fact that such a 2D model gives the same resuh for the 
interaction corrections as the tight binding model on the 
honeycomb lattice with zero spin-orbit coupling. 

In Section III the correction to the self energy of a 
quasiparticlc with momentum p is considered. It is shown 
that the Fermi velocity renormalization in 3D (already 
noted in[lj]) is logarithmic in the UV cutoff A ^ n/a {a 
being the lattice spacing ) very much like in 2D|2l|: 

£d + Aer 
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(4) 



Section IV is devoted to a general derivation and appli- 
cation of the dynamical approach to the electric response 
to the DC and AC electric field. In particular, we obtain 
the slow current decay in a DC field E of the neutral, 
noninteracting 3D Weyl plasma and show that the relax- 
ation is powerwise, see Fig.l. The long time asymptotics 
is oscillating and depends on the microscopic details via 
cutoff a: 
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Here Ny/ is the number of Weyl fermions. This explains 
how the pseudo-Ohmic DC conductivity vanishes in 3D. 
The relaxation dynamics therefore is qualitatively differ- 
ent from 2D, where it is insensitive to the cutoff [2^]. For 
the AC electric field similar slow convergence to the AC 
conductivity occurs, see Fig. 2. Even in the free Weyl 
semimetal one gets, in addition to a finite real (pseudo- 
dissipative) part linear in frequency [lj|, 

an imaginary part, logarithmically divergent as function 
of the UV cutoff: 






(7) 



This is again different from graphenc at zero doping and 
leads to important dielectric properties. Relaxation to 
the asymptotic behavior is faster at higher frequencies. 

The leading interaction correction to the real and imag- 
inary conductivities per Weyl point are subject of Section 
V, see Fig. 3 for comparison with a metal and a semi- 
conductor. The result for the real part expressed via 
the renormalized Fermi velocity defined still linear in fre- 
quency: 
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The normalization frequency is 17 = vp/h. The imagi- 
nary part gets further logarithmically dependent on cut- 
off corrections. 
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despite an apparent " renormalizability" of the model to 
the two loop order at least. This is due to the fact 
that the complex conductivity is not simply related to 
basic Green's functions like in the relativistic theory. 
The physical significance of the results including formu- 
las including that for the complex dielectric constant, see 
Fig. 4, are summarized in the concluding Section VI. We 
speculate about obvious improvements like the random 
phase approximation, renormalization group and a pos- 
sibility of stronger coupling effects like the exciton con- 
densation. 



TIGHT BINDING MODEL WITH DOMINANT 
SPIN-ORBIT INTERACTIONS 

Noninteracting Hamiltonian and linear response 

The noninteracting tight binding model is defined on 
the hypcrcubic lattice n = X)i=i "^i^i with lattice vectors 
a^ of length a by the Hamiltonian: 



^mc = 9 X! rn.icJ,CriCn+a, + ^C. 



(10) 



Here Oi are Pauli matrices, operators d^ create a two 
component spinor a = 1, 2 and Fn^^ is the hopping inte- 
gral that in the presence of an external electromagnetic 
field, described by vector potential Ai, is 



„,j = 7 exp 



/ Ai (n-|-saj,i) 



. ea 



(11) 



where the hopping energy 7 is of order of the band width. 
It is important to derive the current density directly in 
tight binding model: 



J,{r,t) 
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(12) 
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^ / 5^ (r - n - ssii) r„,i4crjc„+ai + he. 



It defines the UV regularization of the current operator 
that obeys the Ward identities. In linear response the 
current density operator is expanded up to the first order 



A as J = Jr 



^f(r) 



Jf(r,i) 



ev 

Y 



^ T„ + he, (13) 
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*|^A(r,i)^T„ + /ic, 

n 

/ 5^(r-n-.sai)c]^crjC„+ai (14) 

Js=0 



where v = "/a/h. Space averages over volume V for a 
homogeneous vector potential A (t) suTiplify: 

jT = ^ / -^f (r) = I^ E 4'T.c„+a. + he, (15) 

^ n 

Expansion of the minimally coupled Hamiltonian in elec- 
tric field is 



K„ 



K + He^t ; H,,t ^-- fjP-A{t). (16) 

C Jr 



Defining Fourier components by c^ = 

A/""^/^ X) s"'*' "cjji where the number of unit cells 

k 

is A/" = V/a^ , one has: 



A' = 7 E sin (ha) c^a^Ck. 



(17) 



Diagonalization of K is achieved adopting a reinterpre- 
tation of the absence of an electron in the valence band 
as a hole and using units ?i = a = u = 1 by 



c£ = v^ak + uZb\ ^ c^+ = v^*a+ + ^r^-k, (18) 



with spinors Vk and Uk given in Appendix A. Up to an 
additive constant it becomes 



K = J2 s^{a 
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where 



Ek = V fcl 



^2 ^2 

^2 + fca , 



(20) 



and the notation k = sin /c was introduced. In 3D one 
observes 8 Weyl points at which £k = inside the Bril- 
louin zone (BZ). Four are right handed: one in the cen- 
ter (the r point) and three on the faces (A"), while four 
are left handed: three on the edges (M) and one in the 
corner {R). The chirality is determined by the expan- 
sion of the Hamiltonian Eq.(fT7| around a Weyl point W, 
sgn{£ijkQ}Q'jQl), where Q* = Jj^jIw- Similarly, in 2D 
one has two right handed Weyls at F and R and two left 
handed at M. 

The paramagnetic and diamagnetic parts of the cur- 
rent due to an electric field (along, let's say, the z direc- 
tion), A (t) = (0, 0, A (i)) ; E{t)^ -~(^li^ (*)' ^'^^■ 
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(21) 



where the functions Xk ^-nd tk are defined in Appendix 
A. 



Coulomb interaction 

Using the expression for the particle densities in mo- 
mentum space (Einstein summation implied for spins 
only) 



kl 



the Coulomb interaction takes the form 



(22) 
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(23) 
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' ~ ^^zZ^p'^i-k-p'^i'-k'-Kp, 
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(24) 



with Fourier transform of the interaction being Vp = 
%^. It is important to note that the electric charge 
of the effective tight binding model should include the 
contributions to the screening due to the polarization 
constant e caused by degrees of freedom not included 
in the model, so the "bare" charge includes this effect 
e^ = e1i/e. We require charge neutrality that is achieved 
by leaving out all the contributions including Wp=o- This 



prescription is always implied in what follows. It is con- 
venient to normal order V with respect to operators a 
and b diagonalizing the " kinetic" term K via Eq. ((T8| : 



V = V 
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(25) 



where the part V''^ contains i creation operators a+ or 6+ 
and j annihilation operators a or 6, all specified in Ap- 
pendix A. In principle the quadratic pair creation V^'^ 
and the pair annihilation terms V'^'^ a could have ap- 
peared. The fact that they have not, explained in the Ap- 
pendix A, greatly simplifies the calculation and makes it 
competitive (at least to the two-loop order) with the dia- 
grammatic approach. The "time independent" approach 
is however much more transparent, when one realizes that 
excitations can be created in fours rather than in pairs 
in this particular model. 



QUASIPARTICLES AND RENORMALIZATION 
OF THE FERMI VELOCITY 



The energy of an electron above the Fermi level with 
quasimomentum p, |p) ~ a+ |0), in the noninteracting 
model described by the Hamiltonian Eq.(??) is: 



p \K\ p) = (O Op ^ £1 (flj^ai + b^bi) a. 



=£. 



The interaction correction is 

A£p = (p|T/|p) = (0|apl/a+|0) 



(26) 



(27) 



Obviously only contributions with equal numbers of cre- 
ation and annihilation operators, V^^ and V^^ can con- 
tribute. The first term contains (see Appendix A) 
a'^b'^ab, b'^b'^bb and a'^a^aa terms, of which only the last 
one could contribute to the expectation value Eq. ip?)) . yet 
it vanishes: 
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We are left with the simple quadratic part V^^: 

qi 
= o / ^^'p-q.9qp' 



where 



(29) 
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This results in: 
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It is instructive to estimate it using the expansion 
around any of the eight Weyl points with the momenta 
restricted, say, by g < A < 7r/2a, thus omitting some 
contributions far from the Weyl points that innocently 
could be thought to be small, but, as the graphene ex- 
ample taught us, might become perfidious if powerwise 
UV divergencies appear in the intermediate stages of the 
calculation. In the present case divergencies are just log- 
arithmic though and one proceeds by expanding around 
one of the Weyl points. 



.9pq ~ cos 0q cos 9p + slu 9q siu 9p cos I 
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(32) 



where spherical coordinates, 

q =(j (sin6'qCOS(/)q,sin0qSin(/)q,cos0q), were used. 
The sum is transformed into an integral: 
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The integrals over angles can be performed, see Appendix 
B, resulting in (for p =p (1, 0, 0)) 
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where r = q/p. Changing variables in the remaining 
integral establishes the leading linear dependence on p: 
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Comparing with the exact result for the "universal" 
model, Ea.([5T|). we can choose the value of the cutoff 
to be slightly outside of the continuum model applicabil- 
ity range: A = 0.53- . This means that the continuum 
model cannot determine correctly the constant term that 
is model dependent. To summarize, one indeed observes 
just a logarithmic divergence and therefore can renormal- 
ize the Fermi velocity written in physical units in Eq. @ . 
In 2D the corresponding calculation gives a well 
known " running" of the graphene velocity towards higher 
velocities [2l[. Therefore the running of the Fermi veloc- 
ity in Weyl semimetal is very much like in (undopcd) 
graphene. 



DYNAMICS OF ELECTRIC RESPONSE: 

CURRENT DECAY IN DC FIELD AND 

INDUCTIVE RESPONSE IN AC FIELD OF THE 

NONINTERACTING WEYL FERMIONS 

NEUTRAL PLASMA. 

In this section we first calculate the time evolution of 
the current when the external electric field is switched 
on in order to obtain the DC and AC electric response. 
This is done however within a "time independent" for- 
malism to carefully trace the states that contribute to 
the dynamics and investigate the emergence of the steady 
state. 



The electric current evolution wthin linear response 



We use the dynamical approach[22| to the semimetal 
response rather than the more customary diagrammatic 
method not just to investigate dynamics under DC or 
AC field, but also to clarify several fundamental is- 
sues in the following sections concerning the Coulomb 
interaction corrections to the AC conductivity. First 
the nature of the dependence of the physical quanti- 



ties on the UV "cutofi" 1/a (renormalization) is eluci- 
dated, next the time evolution is exploited to demon- 
strate the inductive response of the 3D semimetal even 
without electron-electron interaction and contrast it with 
the purely pseudo-Ohmic response of graphene. 

The paramagnetic contribution in the Heisenberg pic- 
ture is given by the vacuum expectation value (VEV): 



(JP(r)) = {ij„,c{t)\3P{r)\ij,nc{t)) (36) 

= (^(0)|C/-i(i)J^(r)C/„,(0|V'(0)). 

The ground state ip (0) is that of an interacting electron 
system without the external field. Expanding to "linear 
response" in coupling to the external electric field (con- 
sidered homogeneous and oriented along the z direction) , 
Eq.dlll), 
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where Ur it) — e *(^+^)*j one gets 
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(38) 



The diamagnetic current, Eq. p^ is already of the first 
order in electric field. In the next section the full model 
including the many body effects will be considered. In 
this section we neglect the interaction V to solve exactly 
for the time evolution of the current density. 
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(40) 
Similarly the diamagnetic contribution, using Eq. (|21[) , is 



Time dependence of the current generated by a time 
dependent electric field 



In the absence of the Coulomb interaction the VEV 
of the average paramagnetic current density (in the field 
direction z, see Eq. (fT3|) ) takes a form: 
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To continue one has to specify the time dependence of 
the applied electric field. We start in the next subsection 
with the constant electric field and then continue to the 
AC case. 



Decay of current in the DC electric field 



Using definitions of the tight binding Hamiltonian and 
the current density operator, Eqs. (fT9| and (|2T|) . one 
obtains 



In the homogeneous DC electric field described by the 
vector potential A (f) = —cEt, one obtains from Eg. pO)) 
the average current. The integral over ii results in: 
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FIG. 1: Fig.l The response of the free Weyl semimetal to a 
DC electric field, multiplied by t, is shown as a function of 
time. The current's asymptotic value of zero is approached 
powerwise. 



in graDhenej22|. Since the BZ can be taken periodic in 
the quasimomentum component kz , the integral over the 
derivative vanishes. The physical arguments put forward 
in the framework of graphcnc[20j in order to comprehend 
this apply equally well here and are not repeated. The 
remaining non-accelerating part, 
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(44) 
is presented in Fig.l. It decays as given (in physical units) 
in Eq.(l5]), exhibiting the zero conductivity rather than a 
universal finite value as in graphene. This asymptotic 
value of zero is approached therefore powerwise and os- 
cillating. Note the dependence on the ultraviolet cutoff, 
that did not appear in graphene|24l|. 
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One expects that the neutral plasma system that does not 
possess electric charges on the Fermi level will not have 
the "acceleration" terms linear in time t that appear in 
the above equation. Indeed the sum of the first two terms 
vanishes: 
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This is seen as follows. The integral over the BZ can be 
represented as an integral over the full ^j^ derivative like 



0. 



Dissipative and inductive response to the AC 
electric field 



The AC electric field, E cos i^lt), is represented by the 
oscillating vector potential, A{t) — —^ sin (fit). The 
" universal" nearest neighbours tight binding model is rel- 
evant only for frequencies n smaller than the hopping 
integral tM- However, as we noticed in the context of 
graphene [20|, the use of the expansion of the dispersion 
relation near the Weyl points should be done with some 
care. 



Performing the integral over ti in Eq.(|40|). one now 
obtains the current: 



J 



JO it) 



E vn'^ 

k 



e^'kism{nt)-2e^' 1 



kl 



n sin (2ekt) - 2ek sin {Q.t) 1 
n^ - Ael J 



(45) 



r 



As in graphene, it is very difficult to approach frequencies numerically when one subtracts the vanishing accelera- 
of interest Q, « 7/fi., when the conductivity is presented tion term Ea. (|43p . that we have encountered in the DC 
in this form. However it becomes substantially simpler calculation, with t replaced by 1/51: 



jo jt) _ e2 ^ 
E " V^ "^ 

k 



l-k'. 



Vt sin (J7i) — 2ek sin (2ekO 



f72 



44 



(46) 



This is shown (together with the electric field) in Fig. 2 certain relaxation time the response becomes periodic ex- 
for several values of frequency. One observes that beyond hibiting a phase difference between the current (points) 
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FIG. 2: Fig.2 The response of the free Weyl semimetal to an 
AC electric field (green lines) is shown as a function of time 
for two different frequencies (magenta lines). The value of 
vja is typically of order ~ 3 • VS^^Kz. 



and the electric field (solid lines). To obtain the steady 
state value of the complex conductivity we average over 
time with a damping factor r] using 



(T (51) = lim — / ( 



^i^t — rjt 



Jit)- 



(47) 



The time integration in Eq. (|46|) with the current density 
of Eq. (|46|) results in 



aoin) 



V 



■E' 



1 -k: 



1 



44 



(48) 
where 0-|_ = fl + ir]. The exact integral is presented (in 
physical units) in Eq.(j6]) and Eq.(l7]) with nonuniversal 
cutoff related to the lattice spacing via A = 0.73- . 

Now wc proceed to try to use the continuum model 
near the Weyl points. After the subtraction of "anoma- 
lous" acceleration terms most of the contributions for 
fi << 7//1 come from the immediate neighbourhoods of 
the Weyl points [23. Due to symmetries it suffices to con- 
sider one of them. The neighbourhood of the origin will 
be defined by fc < A < 7r/2. Here two differences emanate 
compared to the transport in graphene. In graphene the 
lattice spacing is irrelevant, so that the conductivity can 
be calculated within the continuum Weyl model and the 
AC conductivity is real for all frequencies fl << j/h. In 
the present case an inductive part appears and moreover 
depends logarithmically on the lattice spacing. To sec 
this explicitly, let us calculate the AC conductivity ap- 
proximately by using the Weyl approximation (lineariza- 
tion of the dispersion relation) for the contribution of 
the Weyl points. Within the continuum approximation 
in spherical coordinates one writes the sum in Eq. (|48p as: 



aoin) = - 



(27r) 
32e2fi 



fc=0 
A 



sin'^6' 



O 



4fc2 



^\ 



(49) 



3(2^)2 7,^o4fc2-f^2^ 



where the factor 8 is due to the multiplicity of the Weyl 
points in our "universal" tight binding model. The inte- 
gration over fc finally gives: 



ao(r!) 



e^f] / i 4A2 - 



r!2 



(50) 



Stt 



■log 



4A2 



In the last line small terms of the relative order fi^/A^ 
were omitted. The finite result for the real part and 
regularized Matsubara conductivity were first obtained 
in [14 1 . 

Let us emphasize that the electric properties in 3D 
differ from that in 2D (graphene) in that there appears 
an imaginary part of conductivity that is of the same size 
and sometimes even larger than the pseudo-dissipative 
one. Now we turn to the study of the many-body effects 
in 3D Weyl semimetal. 



THE INTERACTION CORRECTION TO THE AC 
CONDUCTIVITY OF THE WEYL SEMIMETAL 

Leading corrections to the current evolution in a 
time dependent field 

The expressions for the diamagnetic and the param- 
agnetic components of the current density as a linear re- 
sponse to arbitrarily time dependent electric fields are 
given by Eqs.(fT3| and (p8| respectively. The evolution 
operator however should be modified as 



^-^I^K+V)t ^ ^-iKt 



T=0 



(51) 
and the ground state (at initial time or Schrodingcr) as 

|^(0)) = |0) + |V'i); l^i) = -;|W4'°M0) (52) 

with the correction to energy being Ai? = (0 |F| 0) =0, 
since the constant term in the energy will be consistently 
omitted after the Hamiltonian was normal ordered. Only 
creation of four particles at once is possible within the 
"universal" model, as is shown in Appendix B. The pair 
creating part V^"^^^' vanishes. This simplifies the problem 
significantly. For example, the diamagnetic component 
is absent due to this. Indeed the diamagnetic current is 
quadratic, Ea (jl3p . and hence the correction 



A/(t) = (V'i 



iKt-diKtl 



^re"'-\Q)+c.c. = {i'l \f\ 0)+c.c. = 0. 

(53) 

The paramagnetic correction (in homogeneous electric 

field along certain direction described by vector potential 

A (t)) of Eq. ([38|) in the presence of Coulomb interactions 

takes a form: 



fit) 

E{t-h) 



A{ti)E{t^h) + cc; 



(54) 



ti=0 



V(0) 
^(0) 



QiHt-p^-iH{t-tt)-p^-iHtt 



V'(O) 



jPQ-iH{t-ti) -p 



V'(o) 



Expanding in interaction with help of Eqs.([5T|) and (|52p . 
the correction to the expectation value in the interact- 
ing ground state S above in terms of the unperturbed 
vacuum expectation values is: 



AS(i-ti) 



AiO 



t-ti 



-iKit-t 



i-r) /^(l.l) _,_y(2,2)\ 



y(2a) ] p~iKTjp 



p^-iK{t-ti)^pJ_y(4,0) 



re 



K 



1 



y(o.4)_LjPg-iA:(i-ti)^-p 



r 





(55 



The other contributions vanish, again due to absence of 
pair creations at this order. The first term corresponds to 



the fcrmion self energy correction, while the rest describe 
the vertex correction. These are readily calculated: 



J 



Aj-(t) 



e 

V2 



"'ti=0 



(56) 



it - U) cos [2s, it - t,)] Ix.r 9.. + ,^-^!i!il?£all^M + W sin[2ep(t-t,)]-sin[2e,(t-t,)] 



where g^p is defined in Eq. (|30p and 
5qp = -Re 






.9qp = Re 



x^x*p{vq-vpy 



This can be calculated using the same methods as in the 
leading order in the previous section. The DC conduc- 
tivity correction vanishes linearly. Finally we calculate 
the correction to current for the AC electric field. 



Correction to the AC conductivity 



In the homogeneous AC electric field, A (i) = 
— -^ sin iilt), the AC conductivity averaged over time as 
in the leading order, Ea. (|47)) . is: 



J 



aiin) 



2ie^ ^ iJp_q 

V2f7 ^ Ael - nl 
pq q + 


IXql'3qp(44 + ^') 2.g-p£q g+^ {As^Ep + n^) 

^ei~n\ sq + sp ^ei-n\ 



(58) 



r 



where r2+ = il + i?;. Subtracting the DC limit (that order, one obtains a much more converging expression 
vanishes after averaging over time) like in the leading proportional to the frequency: 



ai{n) 






E 






|Xql^9qp(l2e 



r!2 



25, 



^qp 



e\{^e\-^\) £q(£q + £p) 



£q£p (,4ep 



■ ^ei + 4e^ 



"D 



r!2 



10 



(59) 



The numerical evaluation results in (without subleading nates around one of the points (chosen to be the origin 
terms in the imaginary part) 



ox(Sl) 



9^ 



log- 



* 1 2 ■"■ 



r): 

(60) Xk = -ie"^''sin6'k 

5qp = - sin 6'q sin 6lp 



(61) 
(62) 



It is however instructive to obtain this expression within 
the continuum (Weyl) approximation valid within radius 
A around each of the eight Weyl points. 

Near a Weyl point one has the following expansion of complementing Eq. ([32l) for gqp. The sum and the 



the functions appearing in Eq. ((5^ in spherical coordi- 



iq — 0p ) (cos 6'q cos 6'p ± 1 ) + sin 0q sin Sp] 

; Eq.dHg) for gqp. The sui 
Coulomb potential in Ea.([5n)) are written as 



A-Ti-p q sin 6'p sin 0q 



2 {li^f Jpq,p^<i>^epe^ p'^ +q'^ - 2pq [cos 6'p cos ^q + sin 6lp sin 6'q cos (0p - (/)q)] ' 



The integrals over all the angles can be performed, see 
Appendix B, 



J 



(63) 



G 



2ie^n 
'(2^ 
2(1 + 



1 



.0 V - ^l 

\2 



3r 



■log 



(1 



(1 



2G~ (r) Apq + Aq^ + Ap^ 



3' 



G^ 



V + fi2 



n^ 



G+{r) 



2G+ (r) 12^2 _ J72 

r 4g2 + rL\_ 



g^2 



,,ii±zt^,^ni±ry 



(1 



3r 



(64) 



with r = q/p. The integrals give 



<yim 



9^ 



log 



4A2 



z 2 4A2 



(65) 



One observes that the subleading terms cannot be given 
correctly by the Weyl approximation as expected. We 
can use the renormalized velocity from Eq.@ to make 
the real part finite.. 



Renormalization of the perturbative expansion for 
conductivity 

We have calculated the renormalization of the Fermi 
velocity, Eq.(|3]) and the AC conductivity of the clean 



Weyl semimetal within the leading order in Coulomb in- 
teraction. It is very tempting to try to use renormal- 
ization to improve the results of the "bare" perturbation 
theory presented to the " two- loop" order above by renor- 
malizing the parameters of the theory. The expression 
for the AC conductivity both in the leading order and 
for the interaction correction contains dependence on the 
UV cutoff a or A and it is interesting to ask whether phys- 
ically measurable quantities can be rewritten via "renor- 
malized" parameters only or the microscopic details rep- 
resented by the cutoff dependence are indeed unavoidable 
for certain measurable quantities like the AC conductiv- 
ity. 

The possibility of " renormalizability" is expected for 
the 3D Weyl model on the following two grounds. 

1. It has been claimed recently and shown to a very 
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high order expUcitly that the 2D version (graphene) is 
renormahzable[27|. This is surprising due to the break- 
ing of the relativistic invariance (that ensures the renor- 
malizabihty for a 2D model with a 3D electromagnetic 
coupling). The AC conductivity expression in 2D can 
be written via the rcnormalized Fermi velocity, Eq.Q 
replacing the bare one in Eq. (p5|) . The key point here 
is that the leading order conductivity is independent of 
both the velocity and the UV cutoff. 

2. The relativistic version, QED, is renormalizable. 



Breaking of the relativistic invariance by taking just the 
static part of the interaction might not spoil this like in 
2D. 

We perform the renormalization of parameters up to 
the two-loop order. 

Replacing the bare Fermi velocity with the rcnormal- 
ized one from Eq. (J4|) in the real part of the conductivity 
from Eq.®, the two- loop result becomes finite and pro- 
portional to frequency: 



J 



'3mJr 

SnVr 






Stt 






log ^^-5 



(66) 



O a' 



where a = -^ leading to the result (in physical units) 
given in Eqs.®. As in 2D the constant C is positive and 
of order 1 . Note that the coefficient e^ in the definition of 
the current is not directly related by relativistic invari- 
ance to the Coulomb interaction in the present model, 
and thus is not rcnormalized. In the order considered 
there is no need to rcnormalize the static Coulomb cou- 
pling. The imaginary part remains therefore " divergent" , 
namely logarithmically dependent on microscopic details. 



DISCUSSION, EXPERIMENTAL 
CONSEQUENCES 

Before discussing the applicability of the results to ma- 
terials proposed as realizations of the Weyl semimetal, let 
us summarize the electromagnetic properties of the clean 
3D semimetals at zero temperature that can be extracted 
from the AC conductivity, Eas. (|6l7l8l9p . Extension of our 
formulas to finite temperature is trivial via the Matsub- 
ara substitution $7 — > 51 — ih/ksT. The properties of the 
clean semimetal are expected to be dissimilar from those 
of a band insulator and a metal. In particular, the di- 
electric and optical properties differ markedly from both 
of them. The complex conductivity of the neutral Weyl 
plasma at zero temperature is very different from semi- 
conductors and from metals. Let us contrast it with the 
Lorentz model of a band insulator and the Drude model 
for metals with electron density n and relaxation time 
T. The AC conductivity of the Lorentz model of a band 
insulator (semiconductor), represented in Fig. 3 by the 
red line, is 



m 



UJ^T 



n 



where the central frequency of the band was taken to be 
rather small wq = 3 • 10^'^Hz, the relaxation time r = 
2 ■ 10"^'* sec and the "plasma" frequency ujp = 10^'^Hz. 
The Drude conductivity of a metal, represented by the 
blue line, is obtained from this formula by taking ujq = 
and the values of r = 2 • 10^^'* sec and Up = 2 ■ 10^'^Hz. 
These are compared with the real part and the imaginary 
part. Figs. 3a and 3b respectively, of the Weyl semimetal, 
represented by the magenta line. The number of Weyl 
points is Nw — 8, the UV cutoff A = n/a, a ~ 3A, 
V = lO^m/s and the intrinsic dielectric constant e = 3 
(due degrees of freedom not included in the model) . At 
low frequencies the absorptive part is linear, so that at 
DC the Weyl semimetal is insulating. However it be- 
comes comparable with that of (a poor) metal at THz 
frequencies. The imaginary part is linear as in a metal 
but with the opposite sign (capacitive like in insulator 
rather than inductive). 

The complex dielectric constant of the 3D semi-metal 
is £ = 1 + Ania/n = e' + ie" 






Nwe^ Aw 
2)-Kvh fl 



a Aw 



(68) 



6vh 



1+ a 



2 , f^ ^ 
3^1°^^ + ^ 



47r n + iriujl-n"^)' 



(67) 



where a = e^/evh with the rcnormalized (measured) 
value of the Fermi velocity v. The normalization fre- 
quency is O and C = —^ ~ —0.53. The real part of 
the dielectric constant e' > 1 is like that of an ordinary 
dielectric material with a very weak frequency depen- 
dence, see Fig. 4a, despite the nonzero AC conductiv- 
ity which is linear in frequency. Albeit note the loga- 
rithmic divergence for small frequencies: for example at 
51 = IQQMHz, e' = 16, see Fig. 4a. The imaginary part 
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FIG. 3: Fig. 3a The real part (Fig. 3a) and the imaginary part 
(Fig. 3b) of the AC conductivity of the Weyl semimetal, (ma- 
genta hne) are compared with those of a band insulator (red 
line) and of a metal (blue line). The parameter values are 
given in the text. 



£ '>-< 




FIG. 4: Fig.3b 



of the dielectric constant depends also only weakly on 
frequency and is universal in that it is of order 0.2Nw 
for Fermi velocity v ~ lO^m/s and intrinsic dielectric 
constant e = 3. 

Let us compare the electric and optical properties of 
the pure 3D Weyl semimetal at neutrality point with the 
corresponding ones in 2D. In the pure 2D Weyl semimetal 




FIG. 5: Fig. 4a The real part (Fig. 4a) and the imaginary part 
(Fig. 4b) of the dielectric constant of the Weyl semimetal, (ma- 
genta line) are compared with those of a band insulator (red 
line) and of a metal (blue line). The parameter values are 
given in the text. 



FIG. 6: Fig.4b 



at the neutrality point (suspended undoped graphene is 
considered to be a good realization of this model [8[) the 
real part of the conductivity is finite and frequency in- 
dependent in infinite samples and with no contact work 
function 2^ , and there is no imaginary part in the AC 
conductivity, at least to leading order in interactions. 
The transport therefore is purely pseudo-dissipative. As 
noted before, the situation changes rather dramatically 
in 3D. We expect that problems of separation between 
the interband transitions (between the valence and the 
conduction bands or electron-hole pairs effects) consid- 
ered in the present study and the intraband transitions 
(including the Klein scattering) due to potential barriers 
and mesoscopic effects are less pronounced in 3D com- 
pared to 2D. Let us now discuss the limitations of the 
model and point out some immediate extensions. 

A clean system was assumed, while disorder is expected 
to be present. The effects of disorder, neglecting the 
interactions, were studied in ref.flj]. There might be 
an interplay between the disorder and interaction effects, 
but the basic physics is expected to be unaltered unless 
interaction or disorder are strong. One also can hope 
that, like in graphene, the importance of disorder might 
be reduced compared to expectations based on physics of 
ordinary "nonrelativistic" quasiparticles. The use of an 
approximation like the tight binding model of section II 
or even a continuum low energy Weyl model is justified as 
long as the frequencies considered are much lower than 
the band width. As we argued in section V, the tight 
binding model is rather universal, but the use of the 
effective massless Weyl theory on the condensed matter 
scale has to be dealt with care. 

It is well known in field theory [25J and in graphene 
that massless fermions cause the absence of a perfect 
scale separation between high energies (on atomic scale 
7) and low energies (effective Weyl theory on the con- 
densed matter scale << 7). It was demonstrated in 
the context of graphene[20j that some aspects of the lin- 
ear response physics, including the Coulomb interactions 
corrections [19j, are not dominated by the Weyl points of 
the Brillouin zone at which the effective low energy model 
is valid. For example, large contributions (infinite, when 
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the size of the Brillouiii zone is being considered infinite) 
to the conductivity from the vicinity of the Weyl points 
are cancelled by contributions from the region between 
them. Another famous consequence of this scale nonsep- 
aration is the "species doubling" of lattice fermions[25j, 
which in the context of graphene means that there neces- 
sarily appears a pair of Weyl points of opposite chirality. 
The UV regularization of the effective theory does matter 
and, if one were to use such a model, the only regulariza- 
tion known to date to be consistent with the tight binding 
is the space dimensional regularization developed in ref. 

There are a number of quite straightforward exten- 
sions of the leading order interaction calculation of the 
AC conductivity performed in the present work. Con- 
verting it into the RPA - like approximation is simple, 
but in addition taking into account the finite momen- 
tum transfer is more involved than the analogous calcu- 
lations in graphene[18|. This would allow to study the 
plasmons' effects and even the strong coupling phenom- 
ena like the exciton condensation[28J. Beyond the linear 
response, phenomena like the nonlinear I-V curves due 
to the particle-hole (Schwinger) pair creation and relax- 
ation due to their recombination studied in the context 
of graphene [29j also can be extended to 3D in a straight- 
forward fashion. 

Acknowledgements. We are indebted to T. Maniv, 
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APPENDIX A. THE UNIVERSAL TIGHT 
BINDING MODEL 

Diagonalization of the noninteracting tight binding 
model and expression for the current 

The noninteracting model is diagonalized in 3D with 
the following coefficients of the Bogoliubov transforma- 



tion Eq.dlHl): 



Uk = 



Vk 



W2efc (sk + kA 



2ek (efe + k,) 






^k + kj \ 

kx -\- ifvy J 



(69) 



where k = sinfc and e| = fc^ + fc^ + fc^. Using matrix 
elements 



uIcr3Vk = ; vTa^Vk = -uTa^Uk = — , 

efc £k 



(70) 



one obtains the coefficient of the electric current as 



Xk 



(ky - ik.A coskz 

£k 



tk 



kz cos kz 

2£k 



(71) 



Normal ordering of the interaction 



The normal ordering is quite straightforward. 



V^ 



klk'V 



Vklk'l'0'\ 



(72) 



where the part V"^^ contains i creation operators a+ or 5+ 
and j annihilation operators a or h. The quartic terms 
are: 



J 



O 



,40 



{vi ■ Uk) [vi, ■ Uk')alai;bZf,b 



(0"4) + 



(73) 



031 ^ (013)+ ^ (^,* . ^^) [(^* . ^^,)a+b\a+ak. - {up ■ u^^) a+b\b\,b_i, 
+ {vf, ■ Uk') [(vf ■ Vk) ajal,b\,ak - {u*i ■ Uk)] a'^btkb-ibtk'i 



(74) 



O^^ = - {vi, ■ Vk') {vi ■ Vk) alal.akak' + {u*i ■ Uk) (w;* • Vk') af.btkak'b-i + (u;* • Uk') {v*i ■ Vk) a^btk/O-kb^i' (75) 
- {ui ■ Uk) {u*i, ■ Uk') btkbtk'b-ik' - {vi ■ Uk) {u*i, ■ Vk') afbtkak'b-i' - {u*i ■ Vk) (w;* • Uk') afib^k'^kb-i- 

I 



In principle there are three possible quadratic terms 
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020 ^ (^02^^ ^ ^^_^, (^* . ^^) (^* . ut,,)a+b+_u, - Si-k' (^r • Uk) {vp ■ Uk')a+h+_k, 



O 



11 



(76) 



5k-i' [v*, ■ Vk') {vi ■ Vk) afak- + 5i^k' (ui ■ Uk) {u*i, ■ Uk-) btkb-i' 
- (u;* • Vk) {vi, ■ Uk') {Sk-i'btk'b-i + (5;_fe/a+afc) . 

I 



Relations between the scalar products for unequal mo- Most importantly, summing over momenta, the pair cre- 
menta that appear here ation term vanishes: 



{u*k ■ ui) = {v*k ■ vi)* ; {v*k ■ui) = - {u*k ■ vi)* = ~ (wj* • Uk) 

(77) 
are useful in summation over momenta incorporating 
Eg. (1241) . in particular for two equal momenta 

Vklk'k — -r:Vl-k5l-k'] Vkk'k'l' = T^Vk' -k^l' -k- (78) 
I 



V 



20 



kl 

Y^ vi-k « • Vk) {{ui ■ Uk) - iv*k ■ vi)) btiaf = 0. 

I 



kl 



This leads to numerous simplifications. 



(79) 



APPENDIX B. ANGLE INTEGRALS 

APPEARING IN CALCULATIONS OF SELF 

ENERGY AND CONDUCTIVITY 

Self energy 



The integral over (j)q in Eg. ([33)) for the self energy is, 
using 

I 



^'^ c + dcosA _ 27r 
A=o "^ + fecos A b 



{bc-ad) /\Jd^ -b'^ +d , 

1 /■■^'^ cos 6*^ cos 6'p + sin 6*^ sin 6'p cos ((/)p — 0q) 

27r J^^=o 1 + ?'^ — 2r [cos 6p cos 6q + sin Op sin 6q cos (0p — 4>q) 



(1 + r^ — 2r cos 6'p cos 6'q ) — (2r sin 6'p sin 0g) 



1 



r 



The integral over the azimuth angle for 9p = (our choice leading to Eq. ([M)) . 



(80) 



(81) 
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Corrections to conductivity 

First contribution to the AC conductivity, Eq. ((5^ , in- 
volves the angle integral 



h = 



-j 

Integration over (j)p and A(/) 



sm 



sm fp sm Oq -^ 



\coii On cos { 



sml 



- q^ — 2pq [cos 9p cos ^5 + sin 9p sin 0^ cos ( 



h - <Pq)\ 



using Ea. ((5n|) gives 



(83) 



h = 



27r2 



'^ou,y 



1 + r" 



\/(l+^ 



-1 



2r cos 9p cos 0,) — (2r sin 0p sin 6'g) 



(84) 



where function G^ is given in Ea.(|Ml). The two other 
angle integrals are done in the same way. 
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